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Abstract 

An ensemble of resonance atoms is considered, which are doped into a medium 
with well developed polariton effect, when in the spectrum of polariton states 
there is a band gap. If an atom with a resonance frequency inside the polariton 
gap is placed into the medium, the atomic spontaneous emission is suppressed. 
However, a system of resonance atoms inside the polariton gap can radiate when 
their coherent interaction is sufficiently strong. Thus the suppression of spon- 
taneous emission for a single atom can be overcome by a collective of atoms 
radiating coherently. Conditions when such collective effects can appear and 
their dynamics are analysed. 



PACS: 42.50.Ct Quantum statistical description of interaction of light and matter 
42.70.Qs Photonic bandgap materials 
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1 Introduction 



Inhibition of spontaneous emission has first been predicted [1,2] to occur in three- 
dimensional periodic structures, in which, due to periodicity, an electromagnetic band 
gap develops. As a result, the spontaneous emission with a frequency inside the band 
gap is rigorously forbidden. This type of matter, where photon band gap arises because 
of the structural periodicity in real space, has been called photonic band-gap materials 
[3]. 

Band gaps also appear in natural dense media due to photon interactions with 
optical collective excitations, such as phonons, magnons, excitons, and so on [4,5]. 
This kind of a gap is termed the polariton band gap. If a single atom is placed in a 
medium with the polariton band gap, and the atomic resonance frequency lies inside 
the gap, then the atom spontaneous emission can be suppressed [6,7]. 

Physical reasons for the effect of suppressed spontaneous emission are rather sim- 
ilar for both types of materials, the difference being in the nature of scatterers light 
interacts with. In artificial photonic band-gap materials, a suppression of the photon 
density of states over a narrow frequency range results from multiple photon scattering 
by spatially correlated scatterers. In natural dense media, such as dielectrics or semi- 
conductors, a frequency gap for electromagnetic modes develops as a result of photon 
scattering by optical collective excitations. The suppression of spontaneous emission 
can be explained as follows. Let us imagine an excited atom in a medium, with the 
atomic transition frequency within the prohibited gap. The atom tends to become 
deexcited by emitting a photon. However, since the propagation of photons inside the 
gap is prohibited, the emitted photon is scattered back and is again absorbed by the 
atom. This means that the atom cannot become deexcited or, in other words, that 
spontaneous emission is suppressed. Hence, this effect is characterized by the condition 

lim s(t) > -1 

for the average population difference s =< of > that is a statistical average of the 
Pauli operator of. In the described physical picture, the limit t — > oo is, of course, 
conditional simply implying that the suppression of spontaneous emission happens 
during times much longer than the emission time 7" 1 , where 7 is the spontaneous 
emission rate. The complete suppression for infinite times is, certainly, never possible 
since there always exist spontaneous photons with frequencies out of the polariton 
band gap, so that sooner or later an atom becomes deexcited via nonresonant photons. 
However, for sufficiently large polariton gaps, such nonresonant relaxation processes 
are very slow. The characteristic relaxation time of a nonresonant process occurring 
through a frequency gap A can be estimated [8,9] as 7 _1 (1 + A 2 /7 2 ). Accepting the 
values typical of the majority of semiconductors [4,5,10-12], for the polariton gap one 
has A ~ 10 13 s _1 while the emission rate is 7 ~ 10 8 — 10 9 s _1 . Hence the characteristic 
time of relaxation through nonresonant processes is 0.1 — 100 s, which is essentially 
longer than typical emission times of free atoms. In this paper, we do not take into 
account the nonresonant deexcitation of atoms, which is justified for rather long times 
as estimated above. 

Even though spontaneous emission of a singe atom in a medium is suppressed, the 
situation for an ensemble of atoms can be different. When a collective of resonance 
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atoms is doped into the medium with a band gap, an impurity band can be formed 
within the prohibited gap. This happens in photonic band-gap materials [13] as well as 
in media with the polariton gap [14,15]. For the collection of identical resonance atoms 
with spacing much shorter than the transition wavelength, electromagnetic effective in- 
teraction takes place due to photon exchange. For closely spaced atoms with transition 
frequencies in the gap, the virtual photon exchanges are of such energy that the atoms 
do not experience the existence of the gap [16,17]. When an impurity band is formed, 
electromagnetic radiation becomes possible. Sufficiently strong effective interactions 
collectivize the resonance atoms that start radiating coherently. Thus, the suppression 
of spontaneous emission for a single atom can be overcome by a collective of atoms. 
Hence the suppressed light can be liberated, at least partially. But what would be the 
dynamics of such a collective radiation? 

In the case of photonic band-gap materials, spontaneous emission near the edge of a 
photonic band-gap has been considered [18,19] for a point-like model with the radiation 
wavelength much larger than the model linear size, A L. For polariton-gap media, 
the dynamics of collective liberation has not yet been properly studied. Moreover, the 
consideration of point-like models, with A ^> L, in the case of suppressed spontaneous 
emission cannot be accepted as realistic because of the follows. Assume that an excited 
atom, with a transition frequency within the prohibited photonic gap, is placed into a 
bandgap material. A photon emitted by the atom, before being reflected back, travels 
at the distance called the localization length l [oc that is about several wavelengths, 
hoc > A. It is only when the localization length is much shorter than the sample 
size, li oc <C L, there is sense of talking that light can be localized, hence spontaneous 
emission be suppressed. While for a concentrated point-like system, with A ^> L, one 
has even more k oc > A L. That is, an emitted photon will safely quit the sample, 
being never reflected back. Therefore the localization of light and suppressed emission 
principally cannot occur in a point-like system. 

In the present paper, we consider a system of iV resonance atoms doped into a 
polariton-gap medium, in which spontaneous emission of a single atom is suppressed. 
The general case of arbitrary wavelengths is studied, which makes it necessary to take 
into account local-field effects. The dynamics of radiation is carefully analysed. Two 
basic points difference the present consideration from the works of other authors: First, 
a realistic situation is analysed, when the radiation wavelength can be (and must be, 
for the spontaneous emission be suppressed and light be localized) much shorter than 
the sample size. Second, the transition frequency is assumed to lie not at the edge of 
the gap but deeply inside the prohibited bandgap. The system of units is used with 
the Planck constant h = 1. 

2 Evolution Equations 

Aiming at studying the dynamics of more or less realistic system, we start with the 
Hamiltonian 

H — H a + Hf + H af + H m + H m f (1) 
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containing the following parts. The Hamiltonian of two-level resonance atoms 



1 N 

Ha = ~ E ^o(l+<) , (2) 

1=1 

with the transition frequency ojq. The radiation- field Hamiltonian 

H f = ^ / (E 2 + H 2 ) dr, (3) 

with electric field E and magnetic field H = V x A. The vector potential A is assumed 
to satisfy the Coulomb gauge calibration V ■ A = and the commutation relation 

E a (r,t), ,4V,0] = 4iric5 a(3 5(r - r') , 

where c is the light velocity. The atom-field interaction is given by the standard dipole 
Hamiltonian 

N n \ 

Haf=~J2 (-3 l -A t + d t -E 0l ) , (4) 



in which Aj = A(rj,£); the transition-current and the transition-dipole operators, 
respectively, are 

3i = iu (d<7+ - dVf ) , di = da+ + d*a~ , (5) 

where d is the transition dipole and af are the ladder operators; and E j is an external 
field, if any. The Hamiltonian of matter, H m) can be modelled by a collection of 
oscillators [6,7]. And matter-field interactions are given by the Hamiltonian 

1 Na 

H m f = y~] j m i • Aj , (6) 

c i=i 

where j m j is a current produced by matter at the point rj; N being the number of 
oscillators modelling the matter. 

When considering stationary properties of electromagnetic field, one usually passes 
to the momentum-energy representation by expanding operators in Fourier series. How- 
ever this representation is not as convenient for studying dynamical properties of a 
system of atoms, especially when their transition wavelength is much shorter than the 
system size. Therefore we shall employ here the space-time representation whose brief 
description can be found in book [20] and detailed analysis in Refs. [21,22]. 

The Heisenberg equations of motion, with the Hamiltonian (1), yield for the atomic 
variables the equations 



dt 

daf 



= -iLU a i + (£; d • Aj - id ■ E 0i ) of , 



dt 



2 (M • Ai - id ■ E 0t ) a+ - 2 (A; d* • A, + id* ■ E 0l ) a~ , (7) 
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where ko = uj /c. The Heisenberg equations for the electric field and vector potential 
lead to the Maxwell operator equations, combining which with the Coulomb gauge 
calibration, one gets 



c 2 dt 2 



with the density of current 



N N 

J(r, = E - r + E " 'i) • (9) 

i=l i=l 



The formal operator solution to Eq. (8) reads 

A(r,t)=A mc (r,t)+ - / J (r', t - - |r - r'|) , (10) 



where A mc is a solution of the related uniform equation. Substituting the density of 
current (9) into Eq. (10) results in the vector potential 

A = A vac + A rac ( + A. ma t (H) 

consisting of the potential A vac of vacuum fluctuations, potential A ra d produced by 
radiating atoms, and the potential A mat induced by the medium, 

N 1 / r<- 
A rad {ri, t) = E jj * 



A mat (r t , t) = E jmi * - — ) , (12) 

where the summation excludes the terms with j = i, and the notation = |ry|, = 
rj — rj, is used. It is worth emphasizing that the relation (10) is an exact expression 
immediately following from the evolution equation (8) and identical to the latter. No 
far-zone approximation is involved here. This approximation could be done [23] by 
expanding |r — r'| _1 in the integrand of Eq. (10) in powers of r'/r, assuming that 
r = |r| ^> r' = |r '|. But we do not invoke here such an expansion. Substituting 
Eqs. (10) to (12) into the Hamiltonian (4) and the evolution equations (7), it is 
easy to notice that all these equations contain an effective dipole-dipole coupling of 
atoms responsible for near-zone local-field effects [24]. Expression (10) is not a final 
solution for the observable field but it is an operator relation having sense only in the 
combination with the Heisenberg equations (7). The system of Eqs. (7) and (10) 
follows directly from the initial Hamiltonian (1) and is often employed in optics [20]. 
Let us introduce the quantities 

Ui (t) = < a-(t) > , Si (t) = < a*(t) > , (13) 

in which the angle brackets imply the statistical averaging over the atomic degrees of 
freedom only, not touching the vacuum and matter degrees of freedom. For the double 
correlators, we shall employ the mean-field-type decoupling 

<afaf> = <afxaf> (i ? j) . (14) 
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It is worth emphasizing that, since only the atomic degrees of freedom have been 
involved in the averaging, the decoupling (14) does not kill quantum effects related 
to the vacuum and the matter degrees of freedom, which have not yet been averaged 
out. The variables corresponding to these degrees of freedom will enter the equations 
of motion for the quantities (13) as operator variables or they can be modelled by 
stochastic variables [25]. This is why the decoupling (14) may be called the stochastic 
mean-field approximation. 

Another problem concerns the way of treating the retardation appearing in Eqs. 
(12). We cannot use the Markov approximation just ignoring this retardation since 
this would eliminate all local-field effects that are so important for realistic many-atom 
systems with wavelength shorter than the system size [26]. But we may employ the 
Born approximation 



retaining the account of local fields. Recall that considering resonance atoms, one 
always keeps in mind that the intensities of interactions of atoms with all fields are 
assumed to be much smaller than the transition frequency uj . If this were not so, 
then the notion of resonance atoms as such would loose its sense, because strong fields, 
essentially shifting the atomic energy levels, would completely change the classification 
of these levels and the values of transition frequencies. Under the assumption that the 
transition frequency ujq is the largest energy scale, from Eqs. (7) it immediately follows 
that the approximation (15) is completely justified [20-22]. This can also be called the 
quasirelativistic approximation since in the nonrelativistic limit, when c — > oo, Eq. 
(15) becomes an identity. 

The equations for the quantities (13) follow from Eqs. (7) where we introduce the 
notation for the effective field 



As usual, we also need to take into account the level width 71 and the line width 72. 
The introduction of these widths is done in the standard way [23,25,27] keeping in mind 
that the values of 71 and 72, in general, differ from the emission rate 7 of a free atom. 
This is because the longitudinal relaxation rate 71 is assumed to effectively include the 
influence of matter which atoms are dopped in, for instance, the influence of lattice 
vibrations. The transverse relaxation rate 72 differs from 7 because of the same reason 
of incorporating the influence of matter, for example, by taking account of dynamic 
and, partially, of inhomogeneous broadening. Thus, 71 and 72 are to be treated as 
phenomenological parameters introduced in the commonly accepted way [23,25,27]. In 
this way, we obtain the equations for the average transition amplitude 




(15) 



fi(t) = k < d • A rad (ri, t) > -id ■ E 0i + , 



(16) 



in which the last term is the random field 



£i(t) = k d ■ [A mc (r i? t) + A mai (r;,t)] 



(17) 



dui 
~dt 



(iu + 72) Ui + fis. 



(18) 
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and the average population difference 



dsj 
dt 



2K7« + />») -7i(s»-so) • 



(19) 



The last term in Eq. (19) allows for the suppression of spontaneous emission in the 
case of a single atom. Really, to return to the latter case, we have to put ft = 0. Then 
we get the equation for a single atom, with the solution Si(t) = s , as is discussed in the 
Introduction. But, in the case of an ensemble of atoms, the presence of nonlinear terms 
related to the effective field (16) can make the solutions to the evolution equations quite 
nontrivial. These equations are to be complimented by an equation either for u* or for 
\ui\ 2 = u*Ui. For the latter equation, we find 



Equations (18) to (20) are the main evolution equations we shall consider. These are 
the stochastic differential equations since they contain the random field (17) including 
the vacuum and matter degrees of freedom. Electromagnetic fluctuations coupled with 
collective excitations of a medium are called polaritons [4,5]. Therefore Eq. (17) 
represents the polariton filed acting on atoms. 

Strictly speaking, one should distinguish two different physical cases, depending on 
the relation between the longitudinal relaxation time T\ = 7f 1 and the characteristic 
time of nonresonant processes T n = 7 _1 (1 + A 2 /7 2 ) discussed in the Introduction. As 
is explained in the latter, it makes sense of talking about the suppressed spontaneous 
emission of an atom only for times t <C T n . This suppression can be either dynamic, 
when Ti <C T n , or static, when T\ ~ T n . In the first case, one has to retain the 
parameter 7x in the evolution equations, while in the second case, one has to set 
7i = 0. For generality, we shall consider both these situations starting with the more 
general case of finite 71, from which it is easy to pass to the particular case of 71 = 0. 

3 Scale Separation 

In what follows, we consider the case when an external field is used only for exciting 
atoms at the initial time, but after this the field is switched off, E j = 0. It is worth 
saying a few words about the possibility of exciting resonance atoms whose transition 
frequency lies in the polariton band gap. There can be several ways of doing this. The 
first way is by exciting the atoms with a short pulse of a resonance external field. The 
amplitude of the latter is to be sufficiently large since a weak electromagnetic wave with 
a frequency within the polariton gap cannot propagate through the medium because 
of total reflection. However, when the incident field is strong enough, a monocromatic 
wave can penetrate into the medium with a polariton gap due to nonlinear effects 
[28,29]. Another possibility for exciting the atoms in a polariton gap could be through 
a third level with the related transition frequency lying outside the gap. It could also 
be possible to get the population inversion at the desired level through a set of levels 
outside the gap, by using nonresonance pumping. 



d\ui 



272M 2 + Si(u*fi + f*Ui) . 



(20) 



dt 
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It is convenient to define the arithmetic averages 



i JV 1 N 

«(*) =mT, Mt) , s(t) = - £ *i(t) (21) 
JV 1=1 JV 1=1 



of the variables (13) and also the average polariton field 

1 N 

= T7 E 6(0 • (22) 
iV i=i 

Let us introduce the collective frequency Q and the collective width T, respectively, 

Q = uj + l2 g's, r = l2 (l-gs), (23) 

where the effective atomic coupling parameters are 

37 A sin kpTjj , _ 37 A cos hpr^ 

9 ~ 4 l2 N konj ' 9 ~ 4 l2 N k r l3 ' 1 j 



and the notation 



4 o 2 U!q 

7 = - k^d , k = — , d = d e d 
3 c 



d\u\ 



is employed. Then from Eqs. (18) to (20) in the mean-field approximation, we obtain 

du 

_ = _ + r) u + si + l2 (g + ig')su*e 2 d , (25) 
ds 

— = -A l2 g\u\ 2 - 2{uX + Cu) - li(s - ao)- 
- 2 72 [(9 + ig'){u*e d f + (g- ig'){e* d uY] , (26) 

-2r| M | 2 + s(u*£ + Cu) + l2 s [(g + tg')(u*e d ) 2 + (g - tg')(e d u) 2 } . (27) 

Equations (25) to (27) form a system of nonlinear stochastic differential equations. 
It is worth stressing that the rotating wave approximation (RWA), well-known in res- 
onance optics [20,23,25,27], cannot be blindly applied for simplifying this system of 
equations. Really, the standard usage of RWA is as follows. One assumes that the 
polarization variable u behaves approximately as exp(— iVtt), that is, rotates with the 
frequency Q, while the population difference s does not rotate. With this assumption, 
one neglects in the right-hand sides of the evolution equations the terms that behave 
qualitatively differently to the assumed behaviour of the left-hand sides. Thus, in Eq. 
(25) one should drop the counter-rotating term containing u* as well as the non-rotating 
term s£. Similarly, in Eqs. (26) and (27) for non-rotating quantities, one should neglect 
the rotating terms with u 2 and (u*) 2 as well as with w£* and But, as is evident, 
such a procedure would completely loose information of polariton degrees of freedom £. 
Hence, RWA is not applicable to such stochastic differential equations. Although the 
RWA reasoning is not justified for these equations, one may try to simplify the latter 
by formally omitting only those terms that would be dropped if the polariton field £ 



8 



were absent, but leaving untouched the terms containing £. Then one would come to 
the system of equations 

du ,.„ „, 

-(iQ + T) u + s£ , 



dt 
ds 

— = -^ 2 g\u\ 2 - 2{u*i + ?u) - 7 i(s - so) , 



= -2T\u\ 2 + s(u^ + Cu) . 

dt 

This is yet a system of nonlinear, because of Eqs. (23), stochastic differential equations. 

To solve this system of equations (25) to (27), we shall use the scale separation 
approach [30-32] which is a variant of the Krylov-Bogolubov averaging method [33] 
generalized to the case of stochastic differential equations. To this end, we need, first 
of all, to classify the functional variables of the system (25) to (27) onto fast and slow. 
This can be done by accepting the usual inequalities 



^ « 1 ^ « 1 



< 1 (28) 



and assuming that the average energy of atomic interactions with the polariton field 
is also much less than Q. Then from Eqs. (25) to (27) it immediately follows that 
the variable u is to be classified as fast as compared to the slow variables s and \u\ 2 . 
The next step is to solve Eq. (25) for the fast variable treating the slow variables as 
quasi-invariants, which yields 



u(t) 



rt 

-(in+r)t 



u + s f e^ +r »' at') dt' 
Jo 



e -(W)t ? ( 29 ) 



where u = u(0). This solution is to be substituted into the right-hand sides of the 
equations (26) and (27) for the slow variables, with averaging these right-hand sides 
over polariton degrees of freedom and over explicitly entering time. In accomplishing 
this procedure, we define the effective atom-polariton coupling 



Re 

a = — 

ST r 



Inn - £ < «*(*)£(*) > dt , (30) 



where Re means the real part and <C . . . ^> implies the averaging over polariton degrees 
of freedom. With the form (29) and the condition £(£) 3>= 0, this gives 



Re ,. 1 

a= T 



lim - T dt [ e-^ +r ^-^ < emit') > dt' . 

By assumption, the atom-polariton coupling is weak, 

M<1. (31) 

In order to understand the structure of the coupling (30), let us model the polariton 
field by an ensemble of iV oscillators, so that 

at) =4rE^ i h ^ kt + b l eluJkt ) > ( 32 ) 
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where b k and b\ are Bose operators with the properties 

< b{b k > > = 4fc' n k , < 6 fe 6 fc / > = , 

and the effective width 7 fc is of order 7 outside the polariton gap but is zero inside this 
gap, 

7 fc = , Ui < LO k < LU 2 , Ap = L0 2 - Ui . 

For the polariton field (32), the coupling (30) becomes 



d = — |7a-| 2 
iV v 



n k l + n k 



K-^) 2 + r 2 (uj k + ny + T 2 



(33) 



We assume that the collective atomic frequency Q lies deeply inside the polariton gap 
of width A p . If this width is much smaller than the frequencies at the edges of the gap, 
A p «(!)!, then Eq. (33) can be approximated by the form 



47 

a ~ 



A 2 p + 4r 2 



After substituting Eq. (29) into the right-hand sides of Eqs. (26) and (27), with 
accomplishing the described averaging and introducing the new function 

w = < \u\ 2 > -a s 2 , (34) 

we obtain the equations 

ds , . dw . . . . 

_ = _4^7 2 w-7i(s-s ) , — = -272(1 - gs)w (35) 

for the slow variables. 

Let us note that the system of three equations (25) to (27) has been derived from 
the set of 3N equations (18) to (20), where % — 1, 2, . . . , N, by employing a mean-field 
type approximation for the averages (21). As is obvious, it is impossible to deal with 
a system of 3N equations for realistically large iV — > 00, because of which one always 
has to invoke some approximation permitting one to reduce the system of untreatably 
large number of equations to a treatable finite dimensional dynamical system. The 
most common such a way of reduction is by using the uniform approximation which 
assumes that the variables Ui and s« do not depend on the index i. This approximation 
is known [20,25] to give a good description when there are no electromagnetic spatial 
structures and the considered sample is sufficiently large for boundary effects being 
neglected. The uniform approximation reduces Eqs. (18) to (20) to the same system 
of equations (25) to (27). Therefore, the usage of the averages (21), together with the 
mean-field approximation, is mathematically identical to the uniform approximation. 
However, from the physical point of view, dealing with the averages is more preferable 
since an averaging description serves as a reasonable first approximation even when 
the sample is nonuniform [20,25]. 
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4 Liberation of Light 



Suppose that at the initial time the atoms are excited so that their average population 
difference s = s(0) > 0. For a single atom, when g — 0, from Eqs. (35) one gets 
s(t) = s , that is, the emission is suppressed. For a collective of atoms, when g ^ 0, 
the system of equations (35) becomes nonlinear, and its solutions essentially depend 
on the value of the atom-atom coupling parameter g, defined in Eq. (24). 

According to the notion of suppressed emission, discussed in the Introduction, we 
need, first, to consider the stationary solutions of the evolution equations. For Eqs. 
(35), there are two stationary solutions. One of them, 

s{ = s , wl = , (36) 

tells that emission is suppressed. And another one, 

S .= I, ro . = Ii(f*^i), (37) 
9 ^l29 2 

shows that light can, at least partially, be liberated. The asymptotic stability of the 
stationary solutions (36) and (37) can be studied involving the Lyapunov analysis. 
This is done by calculating the Jacobian matrix associated with Eqs. (35) and finding 
its eigenvalues that are 



A ± = - Ihi + 72(1 - gs)] =f I {hi- 2 72 (1 - gs)f - Z2^w] 



1/2 



These eigenvalues, evaluated at the corresponding fixed points (36) and (37), define 
the characteristic exponents 



= -7i , \ = -2 72 (1 - gso) , 

1/2' 



1 + 8^(1-^0) 
71 



(38) 



The real parts of the expressions in Eq. (38) are the Lyapunov exponents whose signs 
characterize the stability of the corresponding fixed points. 

The Lyapunov analysis shows that if gso < 1, the stationary solution (36) is a stable 
node, while the fixed point (37) is a saddle point. When gso = 1, both fixed points, 
(36) and (37), merge together becoming neutral. In the interval 

K gs < 1 + , 

8 7 2 

solution (36) is a saddle point, while the fixed point (37) is a stable node. When 

7i 

gs > 1 + — , 

8 7 2 

the fixed point (36) remains a saddle point, and solution (37) becomes a stable focus 
with the characteristic exponents 

\± _ ^ 

A 2 — 2 ~F tUJ OSC ) 
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where the oscillation frequency is 



= 1± 

^osc — 2 



- (gs - i) - i 

This frequency defines the asymptotic oscillation period 



nl/2 



2tt _ 4tt 
w osc ^8(^0 - 1)7172 - 7i 

In the case when 71 72, s ~ 1, an d <? ^> 1, the latter simplifies to 




where i\T x = 1 and 72 T 2 = 1, so that T x and T 2 are the longitudinal and transverse 
relaxation times, respectively. 

In this way, if the coherent interactions between atoms are weak, and the atom- 
atom coupling parameter g is small, so that gs < 1, then the system tends to the 
stationary solution (36), that is, there is no liberation of light. When the atomic 
interactions become stronger, so that gs > 1, the system tends to another stable 
stationary solution (37), and a partial liberation of light occurs since 

s* 2 = - < s (gs > 1) . 
9 

The portion of excitation that remains in the atomic ensemble decreases, with increas- 
ing g, as S2 — > 0. The qualitative change of the asymptotic behaviour of solutions to 
Eqs. (35) happens when gso = 1. This equality defines the bifurcation point for Eqs. 
(35), when the dynamical system is structurally unstable [34]. This bifurcation point 
separates the regions where emission remains suppressed (gso < 1) and where light 
becomes partially liberated (gso > 1). For sufficiently strong atomic interactions, such 
that gs > 1 +71/872, the liberation of light occurs by means of a series of pulses that, 
at asymptotically large times, are separated one from another by the period (39). 

To describe the coherent pulse, occurring when ^ = 0, we may consider Eqs. (35) 
omitting there the relaxation term with 71. Then this system of nonlinear equations 
can be solved exactly resulting in the solution 



s = - ^- tanh f LJo) + i , w = -*4 sech 2 (tz*) , (40) 



912 ' " V T J ' g 1 4g 2 >y'£ ~ V r 

where the radiation width 70 is given by the relation 

7 2 = r 2 + 4/ 7 2 2 (l«o| 2 -«o^) (41) 

in which u Q = w(0), a = cc(0), and 

r = 72(1 - ffs ) , 7o = — , 

To 
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and where the delay time is 



to = y In 



70 - r c 



70 + r 

Introducing the critical atom-polariton coupling 



(42) 



(g*o - If 



+ 



PO 



-.2 ' 



the radiation width can be written as 



7o = 2g\s \j 2 Va c ~ «o • 



(43) 



(44) 



The value (43) is termed critical since the coupling a cannot exceed a c for the solution 
(40) to remain finite. The restriction «o < « c specifies condition (31) assumed for the 
validity of the averaging method. This restriction is not as severe as far as « C 1 
while, for large g, a c > 1/4. 

At the time t = t , the solutions (40) and (34) yield 

s(t ) = - , w(t ) = sl(a c - a ) , < \u{t )\ 2 > = sl(a c - a ) + ^- . (45) 
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Then the system of atoms achieves the maximum of coherence. 

Consider more in detail the case when the system of atoms is initially completely 
inverted, sq — 1, there is no triggering pulse, uq — 0, and the atom-atom coupling is 
strong, g > 1. Then the critical value (43) becomes 



9-1 

> 2 9 , 



(46) 



If a <C a c , the radiation width (44) gives 

lo = (9- 1)72 
The corresponding radiation time is 



To 



T = 



9-1 

And for the delay time (42), we have 



1 + 



2ff 2 c*o 



(5 -I) 2 



(47) 



(48) 



to 



To 



In 



When the atomic coupling is strong, then 

1 



(9 - I) 2 ~ 9 2 a 



9 2 olq 



(49) 



a c ~ - (uo = 0, c/ > 1) 



(50) 
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The radiation time (48) for g >> 1 becomes 

r ^ — (1 + 2a ) . (51) 

Hence, r can be much smaller than T 2 . Since g ~ AT, we have To ~ A^ -1 , which is 
typical of superradiance [20,25]. Under condition (50), the delay time (49) is 

T 2 

*o — 75— I ln «o| • (52) 

Note that if ao — > 0, then to ^ 00, which means that emission would be suppressed for 
very long time. The superradiant burst develops at the delay time to and lasts during 
the radiation time r . 

In the case when a c — a <C 1, the radiation time increases, as compared to the 
opposite case when a <C a c , being 

To = 1 1 S = , (53) 

while the delay time shortens becoming 

2g\s \T 2 



(gso - iy 



V a c - «o • (54) 



In this case, the radiation can be coherent if r <C T 2 , which requires g^/a c — a ^> 1. 

The superradiant character of emission is connected with the radiation time being 
inversely proportional to the number of radiators, r ~ A^ -1 . Another characteristic 
describing the level of coherence in the radiating system is the radiation intensity. The 
total intensity of radiation, averaged over fast oscillations, 

/(*) = hnc{t) + hohit) , (55) 

consists of two terms, the intensity of incoherent radiation 

hue = 2 ^07^(1 + s) (56) 
and the intensity of coherent radiation 

hoh = ojqIVsN 2 < \u\ 2 > , (57) 
where ip s is the shape factor [35] given by the integral over spherical angles, 



<P* = -^ J In x e d \ 2 F{k n) dQ(n) , F(k) 



-1 N 2 

— e ik ' ri 



(58) 



For the point-like system, for which k — > 0, we have tp s — > 1, and the intensity 
of coherent radiation is proportional to the number of radiators squared, I co h ~ A^ 2 , 
which is typical of the Dicke model. For a finite-size system, the shape factor essentially 



14 



depends on the relation between the radiation wavelength A and the characteristic sizes 
of the sample [35]. Thus, for a cylindrical sample of radius R and length L, we have 



— <C 1 — 1 

8L > 2-kL ^ L ' L ' 



for pencil-like or disk-like shapes, respectively. Hence for the first case I co h ~ iV 5 / 3 , 
while for the second case I co h ~ iV 4//3 . 

To compare the intensities of coherent and incoherent radiation, it is convenient to 
introduce the coherence coefficient [36] defined as the ratio 

c coh {t) - Hty • ( 59 ) 

For the coherent radiation intensity (57), using the relation (34), we get 

I C oh(t) = uj i<f s N 2 (w + as 2 ) . 
Therefore, the coherence coefficient (59) is 

C coh = 2 Vs N . (60) 

1 + s 

At the moment of the maximal coherence of the superradiant burst, according to Eqs. 
(45), we find 

C coh {t ) = 2g(p s s 2 ) N ac a ° . 

l + g 

If C co h > 1, the radiation is predominantly coherent. For the case when uq = 0, sq = 1, 
and g ^> 1, the latter expression gives 

C coh (t ) ~ i (p s N . 

If the sample has a pencil-like or disk-lake shape, then 



C co h(to) 



ff (f ) 2 pA 3 , pencil 



3 
16?r 



(jj pA 3 , disk , 



where p is the density of resonance atoms. When this density is sufficiently high, and 
A <C R or A <C L, the coherence coefficient can be very large, which would mean that 
the radiation is almost purely coherent. 

The general condition for light to be, at least partially, liberated is 

s(oo) < s • (61) 

From equations (40), we have 

s(oo) = — — + - , w(oo) = , 

912 g 
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which, with the radiation width (44), gives 



s(oo) = — 2|s |v /q; c — «o + - • (62) 

Therefore the liberation condition (61) becomes 

g(s + 2|s |v / «c - «o) > 1 • (63) 

In the case, when uq = 0, and taking into account that a <^ 1, equation (62) reduces 
to 

*(oo) = ( S ,°; / 9 *° < 1 1 (64) 
y J \ -s + 2/g , gs > 1 . K J 

The condition (63) simplifies to 

gs > 1 . (65) 

It is interesting that, although the limit values s(oo) are different for the cases when 
7x is finite or when it is zero, but the liberation condition (65) remains the same. 
Under this condition, a system of atoms can radiate, though spontaneous emission of 
a single atom is suppressed. The case of a single atom can be recovered by setting 
g — > 0. Then, as is obvious, condition (65) can never hold true. It is only when the 
density of doped atoms is sufficiently high, so that the atom-atom coupling g becomes 
sufficiently large, satisfying condition (65), the radiation of atoms is possible being due 
to collective effects. 



5 Discussion 

We considered a system of resonance atoms doped into a medium with a polariton 
band gap. Spontaneous emission of a single atom with the transition frequency inside 
the polariton gap is suppressed. However an ensemble of atoms with their transi- 
tion frequencies in the gap can radiate due to effective coherent interactions between 
the atoms. If this intertaction is sufficiently strong, light is partially liberated. The 
collective liberation of light occurs through one or a series of superradiant pulses. 

The dynamics of light liberation has been analysed employing the scale separation 
approach [30-32], which is a generalization of the averaging method [33] to stochastic 
differential equations. This approach provides an efficient tool for treating complicated 
systems of nonlinear evolution equations, as has earlier been demonstrated for the 
problems of superradiant spin relaxation in nonequilibrium magnets [30-32] and of 
nonlinear dynamics of atoms in magnetic traps [37-39], where nonlinear phenomena 
are of crucial importance. Collective liberation of light is also a principally nonlinear 
phenomenon, with nonlinearity caused by coherent atomic interactions. 

In conclusion, let us give some estimates for the characteristic parameters considered 
in the text. The polariton effect is well developed in many dielectrics and semicon- 
ductors [4,5]. For instance, it is intensively studied in such semiconductors as CuCl, 
CuBr, CdSe, ZnSe, GaAs, GaSb, InAs, AlAs, and SiC [10-12,40]. The polariton gap in 
such materials develops around the frequency 10 14 s _1 , with the gap width A p ~ 10 13 
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s~ 4 . By assumption, the atomic transition frequency is inside the polariton gap, i.e. 
ujq ~ 10 14 s -1 . Hence the radiation wavelength is A ~ 1CT 3 cm. For the line width, one 
may take 72 ~ 10 9 s _1 . 

For the initially inverted atoms, with so = 1, the critical atom-atom coupling 
parameter, above which the collective radiation becomes possible, is g c — 1. For the 
parameters g > g c , the collective liberation of light occurs. If 7 ~ 72, then the atomic 
coupling parameters, defined in Eq. (24), are g ~ pA 3 and g' ~ pA 3 , where p is the 
density of atoms. Consequently, the critical density of atoms, providing g c — 1, is 
p c ~ A~ 3 , which gives p c ~ 10 9 cm -3 . When p > p c , radiation becomes possible 
because of the formation of an impurity band inside the polariton band. The width of 
the impurity band is of the order of the collective line width T ~ 7 2 g ~ 7pA 3 . The latter 
becomes larger than the polariton gap, when the density p > A p /7 2 A 3 , that is, p > 10 13 
cm -3 . For such densities of resonance atoms, the polariton gap can be overlapped by 
the impurity band. At the density p ~ 10 13 cm -3 , the atom- atom coupling is g ~ 10 4 , 
and the collective line width is T ~ 10 13 s~ 4 . 

With the densities p < 10 13 cm -3 , the atom-atom coupling g < 10 4 , while the 
atom-polariton coupling is a ~ 7 2 /Ap, that is a ~ 1CT 8 . The radiation time tq ~ T2/ ' g 
is larger than 1CT 13 s but can be much smaller than T2 ~ 10 -9 s -1 . And the delay time 
to ~ to I In a I is an order longer than the radiation time. In this way, for the density 
of doped atoms p ~ 10 9 — 10 13 cm -3 , the effective atom-atom coupling is g ~ 1 — 10 4 . 
Then the delay time is t ~ 10~ 12 — 1CT 8 s and the radiation time r ~ 10 -13 — 1CT 9 s. 
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